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Secondary nuclear production physics is receiving increased attention given the high-quality measurements of 
the 7-ray emissivity of local interstellar gas between ~ 50 MeV and ~ 40 GeV, obtained with the Large Area 
Telescope on board the Fermi space observatory. More than 90% of the gas-related emissivity above 1 GeV is 
attributed to 7-rays from the decay of neutral pions formed in collisions between cosmic rays and interstellar 
matter, with lepton-induced processes becoming increasingly important below 1 GeV. The elementary kinematics 
of neutral pion production and decay are re-examined in light of two physics questions: does isobaric production 
follow a scaling behavior? and what is the minimum proton kinetic energy needed to make a 7-ray of a certain 
energy formed through intermediate 7r° production? The emissivity spectrum will allow the interstellar cosmic- 
ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input 
to solar modulation models. A method for such an analysis and illustrative results are presented. 



1. Introduction 



The majority of the 'diffuse' 7-ray emission of the 
Milky Way is truly diffuse, with the contribution from 
a superposition of unresolved, low luminosity point 
sources estimated at only the 5-10% level in the Galac- 
tic plane. Almost all of this diffuse emission is the 
consequence of propagating non-thermal hadronic cos- 
mic rays collidin g w ith the nuclei of interstellar par- 
ticles (see e.g. |22[ |23|), and cosmic-ray electrons 
and positrons making bremsstrahlung 7 rays. Be- 
sides making 7-rays a valuable tracer of the interstellar 
matter and radiation, the spectral energy distribution 
(SED) of the diffuse emission encodes the spectrum of 
the interstellar cosmic rays. 



*Also at the Department of Physics, University of Tokyo, 
Tokyo, Japan 



Using data taken in the first six months of Fermi 
science operations, Abdo et al. [l| measured the spec- 
tral energy distribution (SED) of gamma-ray emis- 
sion associated with local neutral atomic hydrogen, 
HI, between 100 MeV and 10 GeV, finding that it 
is proportional to the HI column density and deriv- 
ing an emissivity (emission rate per H atom) spec- 
trum. Many measurements of the HI emissivity in 
the neighborhood of the solar system and the outer 
Galaxy were obtained using Fermi data. More re- 
cently, Casandjian [5] presented a new emissivity spec- 
trum of local atomic hydrogen between Galactic lati- 
tudes 10° < \b\ < 70° for energies - 50 MeV to - 40 
GeV with error bars mainly < 15%. 

It is therefore an optimal time to use the emissivity 
to deduce the interstellar cosmic ray proton spectrum 
with improved accuracy. While our ultimate goal is 
the determination of the interstellar cosmic-ray pro- 
ton and a-particle spectra from the measured diffuse 
Galactic 7-ray emissivity [25j, we first focus here on 
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the topic of near-threshold neutral pion production. 

Due to the dominant abundance of protons in the 
cosmic rays, and plentiful hydrogen in the interstel- 
lar medium, the most important process making the 
diffuse 7-ray glow of the Galaxy is 

p + p^>y + X . (1) 

Here, X is anything else made in the reaction besides 
7 rays, and the implied inclusive description already 
becomes dependent on the environment of the inter- 
action, for example, through 7-ray production by sec- 
ondary leptons. For simplicity of analysis, we sepa- 
rately consider primary and secondary leptonic emis- 
sions, and confine the calculations only to the 7-rays 
formed from the decay of 7r° (8.3 x 10 _1T s mean life) 
produced through nuclear collisions. 

As considered by Neronov et al. pij, the gamma-ray 
spectrum from cosmic-ray proton distributions in the 
form of power laws in kinetic energy exhibits a signifi- 
cant hardening below proton kinetic energy T p = 9^ 
GeV. The importance of the nuclear physics uncer- 
tainties in performing the spectral inversion to charac- 
terize the limits of accuracy of the derived cosmic-ray 
proton spectrum — leaving aside leptonic and metallic- 
ity corrections — may undermine the conclusion about 
a spectral break, especially if the intrinsic spectrum 
is related to a power-law distribution in particle mo- 
mentum, or momentum per unit charge (GV/c) [9[. 

The diffuse galactic 7-ray emissivity provides a 
clearer tracer of hadronic cosmic rays than SNRs, 
which are generally too small to resolve, dim, or in a 
confused region. Recent Fermi observations of SNRs 
IC 443 and W44, and AGILE observations of W44 
[To| . are providing strong evidence for a 7r°-decay fea- 
ture in the spectra of these SNRs [3|- For a general 
description of cosmic-ray origin and propagation, the 
7-ray spectrum made by cosmic rays accelerated at 
their sources should be connected to the spectrum of 
cosmic rays in interstellar space, which is revealed by 
the Galactic diffuse 7-ray emission and, for cosmic-ray 
electrons, by radio synchrotron and inverse Compton 
7-rays. Determining the interstellar cosmic-ray spec- 
trum from the 7-ray emissivity, when related to the 
cosmic-ray spectra in SNRs implied by their SEDs, 
will give information about cosmic-ray escape, and es- 
sential input for testing solar modulation theory. 

Here we gather together a few of the tools needed 
to make an assault on the nuclear physics related to 
the Galactic diffuse 7-ray emissivity from proton in- 
teractions near threshold to a few GeV kinetic en- 
ergy, where resonance production is most important. 
Surprisingly, the radiation physics for cosmic-ray elec- 
trons is known to better than a percent, determined 
by the fine-structure constant, but the pion-decay 7- 
ray production cross section at the GeV-scale is hardly 
known to ±10 - 20% accuracy. This is not likely to im- 
prove until new laboratory physics experiments mea- 



sure these cross sections. In the meantime, we can 
assemble models for comparison. 

Refs. [l9j and list the following near- threshold 
decay products of pp collisions: A(1232), A(1600), 
and N(1440), all of which decay through a 7r° to make 
7-rays. Furthermore, there are baryonic and mesonic 
decays. In the former class are A and A, S° — >• A + 7, 
and E + , with only the E° in a direct, non-7r° 7-ray 
channel. Meson-decay channels other than 7r° are if + , 
if - , Kl, and K$ — )• 27, with only the latter having 
a direct, non-7r° 7-ray production channel. Keeping 
also in mind that direct production occurs in about 
half of the 77 production cross section, which is ~ 10% 
of the total inelastic cross section, a very small cor- 
rection to the symmetry properties of the gamma-ray 
production mechanism is possible. Note the possible 
confusion between direct and 77 production mentioned 
in Ref. 0. 

In this contribution, the kinematical relationships 
of the isobar model are reviewed. For reference, see 
Lindenbaum & Sternheimer [l4| , Stecker [l9| , Dermer 
[6|. We consider whether isobaric production follows 
a scaling like behavior, where the invariant cross sec- 
tion becomes independent of the invariant colli- 
sion energy. We answer what appears to be a simple 
question, namely what the maximum photon energy 
produced through intermediate 7r° production is as a 
function of proton kinetic energy. Then preliminary 
results for the determination of the interstellar cosmic- 
ray spectrum from the Fermi-LAT 7-ray emissivity are 
presented. 



2. Isobar Kinematics 

The proton kinetic energy is T p = m p (j p — 1), 
and the invariant square of the CM energy is s = 
2m p (j p + 1) = 2m p (T p + 2m p ). In terms of the center- 
of-momentum (CM) Lorentz factor 1 

f^. < 2 > 

At threshold, y/s = (2m p + m^o) . Therefore the pion- 
production threshold is 

Tl min = 2m v o{\+ 7 ^-) = 279.6 MeV , (3) 

3.6% larger than 2m^o. The mass of the neutral pi- 
meson is m^o = m n = 134.96 = 135 MeV (m p = 
938.27 MeV). 



1 The (3 factor is related to the Lorentz factor 7 as usual, (3 = 
Y^l — I/7 2 , with corresponding subscripts and superscripts. 
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Figure 1: Kinematics of the elementary p + p — >• 7r° — >• 2j in the Laboratory System (a), in the Center of Momentum 
System (b), and in the CM following the collision (c). Note assumption that outgoing baryons retain same direction as 
the incident cosmic ray. 



In the aftermath of an inelastic nuclear collision, 
the invariant energy remains constant. A two-body 
collision leading to two bodies in the final state is 
a particularly well-defined event, with the final state 
assumed to consist of a baryon in a resonantly ex- 
cited state, which decays into a p and 7r°, after which 
the 7T° — >> 27 decay is essentially instantaneous. Thus 
= m p j*j + rriA^A m pl + m A7A- See Figure HJ 

From momentum conservation, rripfij = t?7a/#a7a, 
implying 7 = yjl + (m^/m p ) 2 (^\ - 1), so that 



7A 7a 



777/ 



mi 



2 v /imA 



(4) 



When the A decays into a proton and pion, Equation 
dH) is still valid in the A rest frame, where yfs = ttia = 



171/ 



mz 



7tt 



and, by symmetry, 



7 P 



2mA^7r 



771/ 



2m/\m p 



(5) 



(6) 



In the "Laboratory System" (LS), the Lorentz fac- 
tors of the forward- and backward-moving isobars are 



7a =7*71(1 ±F /?a)- 



(7) 



The energy distribution of pions in the LS is given by 
noting that the Lorentz invariant 



E 

d^> 



1 



p 2 dpdVt pdEdQ, ' 



(8) 



because pdp = EdE. Thus — — — is also invariant, 

P7 djd/jL 

and 

dN f 1 dN If 1 ,, h , fQ . 
In terms of the 7r° spectrum, 

-<»>■ <io> 

7tt,o = 7^7a(1 ^ /^tt^a)' which depends on the isobar 
mass 777a. Solving the 5 function gives 

dN = { H[^-i^{l- + 

<h* 1 4^; 7 ;/3X7a 



g[ 7 , ; 7^1 7^(1 +/M)] 



} , (11) 



where a, 6) = 1 if a < x < 6, and H (x\ a, b) = 
otherwise. 



2.1 . Cross Section for p + p -» 7t° + X 
Production 

The cross section for 7r° production from p + p in- 
teractions can be written as 



d(J pp ^(T p ) da^Tp) dN^Tp) 



dTir 



dT^ 



dT^ 



(12) 
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Figure 2: Inclusive cross section for 7r° production, from Q. See also [20 



The term {^a^{T p )) is the inclusive cross section for 
the production of 7r° (Figure [2]), irrespective of the re- 
maining content of the secondary beam, and the pro- 
duction cross section is normalized to unity, so 



Jo 



' dNAT E l_ 



By contrast, exclusive cross sections are cross sections 
for specific decay channels. 

The inclusive cross section of particle i in reaction Y 
is the product of the inelastic cross section for reaction 
Y and the multiplicity Q of particle i. Consequently 
the inclusive cross section is much larger than the in- 
elastic cross section at energies well above threshold 
when the multiplicity is high. The inclusive cross sec- 
tion for the production of 7r° in p-p collisions can be 
fit within the uncertainty of the cross section measure- 
ments by the function 



GVxCmb) 321np p 



48.5 



59.5 



(13) 



for proton momentum p p = m p ^/(l + T p /m p ) 2 — 1 in 
the range 8 GeV/c < p p < 1000 GeV/c [7]. See [3 
for 7-ray production cross sections at LHC energies. 



2.2. 7 rays from 7r° Decay 

In its rest frame, the pion decays into two 7 rays 
with energy ra^/2, but in the LS, the 7r°-decay 7 rays 
are radiated with every allowable energy between a 
kinematic minimum and maximum energy defined by 
setting ji = ±1 in the relation e' = m^/2 = 7^(1 — 
/Stt/x). If the 7T° decays isotropically in its own rest 



frame, then the 7-ray decay spectrum in the proper 
frame of the 7r° is 



dN _ 8(e' - rriir/2) 



de'dSl' 



47T 



(14) 



The factor of two arises because two photons are pro- 
duced per interaction. For a 7r° produced with Lorentz 
factor 7^, the transformation properties of N(e, fi) im- 
ply from Equation (|9j) that 



dN 



5(e' - mj2) 



dedn 2^(1-/3^) ' 

so that 
dN 2 

de (3^171^ 



(15) 



H[e; ^77^(1-/^), ^77^(1 + /^)] . 

(16) 

The spectral number emissivity for tt° production 
from cosmic ray protons colliding with target protons 
is 

do- p H^7r (Tp) 



f'OG 

i P H->Tr°(T v ) = 4im p / dT p jp(T p ,Q p ) 
Jo 



(17) 

wherer j(T p ) is the cosmic-ray proton flux, described 
in more detail below. The 7-ray emissivity from 7r° 
decay is 

n (e)= — r dT h P H ^°^ 
«°^ e > m OT J T „ in - v /r„(T„ + 2m ff ) ' 

( 18 ) 

where T™ m (e) = m 7r [7^' m (e) — 1], and 7™ m (e) is given 

by 

7r n ( e ) = J[^T + ^^], (19) 



2 L (m n /2) 
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that is, the solution is symmetrical to the inversion 
x — >• 1/x, where x = e/(m 7r /2). This is the famous 
bilateral symmetry of the pion bump around 67.5 MeV 
in a photon spectrum [19j , which is seen as a hardening 
in a vF v spectrum below several GeV. It is obvious 
from its appearance but requires a moment to see that 
the two-parameter log-parabola function 



da(T p ) 
de 



Kx~ blogx 



(20) 



preserves the x — >• x~ l symmetry, and might provide 
a useful characterization of the p + p — » 7r° 27 
spectrum with, in general, K = K(T p ) and b = b(T p ). 

2.3. Breit-Wigner Distribution 

The resonance mass distribution is assumed to be 
described by the Breit-Wigner function, 



B(m A ) = w r (T p ) 



7r[(m A -m° A ) 2 -f T 2 ] 



(21) 



where T is the resonance width. The normalization 
factor 



w r (T p ) = 7r[arctan 



171/ 



arctan 



rriv 



171/ 



(22) 



2.4. Production Kinematics 



We calculate the production spectrum of secondary 
7r from cosmic-ray protons with flux j p (T p ,Clp), in 
units of cosmic ray (CR) p/(cm 2 s sr GeV), colliding 
with protons at rest. Hence 



poo 

47rn p / dT p j p (T p ,tt p ) 
Jo 



da(T p ) 
dT n 



(23) 



The spectral emissivity is given, in units of 
ph/(cm 3 s GeV H atom), by 



s(e) 



f'OO 7771 pOO 

= 8tt —a- dT p j(T p ,n p ) 



da{T p ) 



(24) 

The minimum proton kinetic energy, T™ m (e), re- 
quired to produce a photon of energy e made through 
intermediate 7r° production turns out not to be ex- 
pressible in simple analytic form, and we calculate it 
numerically below. 
For the isobar model, 



,rr = -j — / dm A B(m A ) -j— , 

dT n Awm n J mv+mn dj n 



(25) 



with dN/dj-n given by Equation ([TT]) . By comparison, 
at high energies, scaling representations are needed to 
describe secondary tt° production; e.g., the Stephens 
& Badhwar [2l| invariant cross section is given by 



E ■ 



d 3 



0~SB 



dp 3 



Af(E p )(l - x)" exp[-Bp ± /(l + 4m 2 p /s] , 

(26) 



where E p is total proton energy, 



4m 



/ )2x(1 + ^ } ' 



f(E P ) = (1 



and q = [c x + p±(c 2 + c 3 p±)]/yjl + 4m 2 /s, with 

A = 140.0,5 5.43, ci = 6.1, c 2 = -3.3, and 
c 3 = 0.6. In Equation (j26j) . £ = y^* 2 + x* 2 , 
where ^ = p£ cos 9/pl iTnax , p± = P*sin<9, 6> is 
the polar angle in the LS, p* 



4, 



£;* = 7*(£; 7r - f3* Pit cos Omax)- The quantities = 



V^bi + m 2 )/s and p£ 



E i 2 max ~ m li where 



(s + m 2 — 4m 2 ) / 2^/i. See Dermer & Menon 



Stephens & Badhwar [21] for more details. 



2.5. Form of the Interstellar Cosmic-Ray 
Spectrum 

Now consider 7r° production from interactions of 
cosmic rays with Galactic gas. The dominant con- 
stituent of the interstellar medium is H, and we con- 
sider only the dominant proton component of the 
galactic cosmic rays. A suggested form for the de- 
modulated cosmic-ray proton spectrum observed in 
the Solar neighborhood is 



j* T *(T p ,n p ) = 



2.2 

£^75 



(28) 



[6|, namely a power-law in total energy E p = T p + 
rripC 2 . In the same units, Dermer [9] suggests that a 
more physical model is a power law in momentum, 



f p lpp (T p ,n p ) 



2.2 



2.85 



(29) 



In Figure [3j we use Equation ([29]) to make cal- 
culations of the emissivity spectrum for the hybrid 
model treated in Ref. [6|, as well as for the iso- 
baric and scaling models separately. Note that the 
Stephens-Badhwar scaling function, even when used 
to model the entire production spectrum, is well be- 
haved, whereas the isobar model overproduces sec- 
ondary production if extended to high energies. 



3. Discussion 

We motivated this study by two questions. The 
first is: does the isobar model follow a scaling rep- 



eConf C121028 



4 Fermi Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



E 

u 

> 



-l-J 

'l/> 

'£ 

LU 

X 




10 



100 



1000 
E(MeV) 



10" 



10 3 



Figure 3: Calculations of p -\- p — » 7r° — »• 2j production, compared with the preliminary Fermi-LAT emissivity derived 
by Casandjian 5], which is shown by the red data points with error bars. The proton flux is given by a broken 
power-law spectrum described by the power-law momentum spectrum of Equation (|29)) with index s = 2.85 below 200 
GeV, which is connected to a power-law spectrum with s — 2.65 above 200 GeV to reproduce the PAMELA hardening 
0]. The model emissivities are multiplied by a factor of 1.84 [3] to account for the composition of ions heavier than H 
in the cosmic rays and ambient medium. The heavy solid curve uses the hybrid model of Dermer [(J, with separate 
isobar and scaling components shown by the long-dashed light purple and short-dashed light blue curves, respectively. 
The blue dotted curve is the same calculation using the Kamae et al. [13] model. An all-isobar calculation and an an 
all-scaling calculation (using the Stephens & Badhwar [2l|] scaling function) are shown for comparison by the 
light-dotted orange and dot-dashed green curves, respectively. The all-isobar component is shown only to illustrate its 
high-energy behavior; the isobar cross section declines oc 1/y/s at high energies. Bremsstrahlung is not included, which 
can also make a significant contribution, especially at low energies. 



resentation? Most simply, the Feynman scaling hy- 
pothesis that Eda/d 3 p becomes independent of s at 
large s ^> m^. By inspecting Equation (|26]h one can 
see that while the invariant cross section becomes s- 
independent in the limit s ^> m^, this s-independent 
behavior is not followed by the isobar model. Figure [3] 
shows that unusual structure appears in the extreme 
situation where every event follows A(1232) resonance 
formation and decay. As can also be seen, structure 
in the production spectrum is very sensitive to the 
relative isobar ic and scaling contributions. 

The second question is, what is the maximum 
photon energy emerging from secondary production 
through intermediate 7r° formation? At threshold, 
Equation (j3]), the maximum photon energy e max is 
clearly At high energies, e max « T p . We 

have derived the equations to answer this question 
precisely. The maximum photon energy 

ir, max ) , (30) 

where, using Equation (fTTjh 

7*,™* = 7^(1 + (31) 

Substituting Equations ((2]) and ((4]) into Equation (O, 
and using this result and Equation ((5]) in Equation 



(|3T|) gives the result in terms of ttia, where the allowed 
isobar mass range 



m p + m Jlx < rriA < \fs — rn p , (32) 



and depends on proton kinetic energy T p . The nu- 
merically calculated value of e max /T p vs. T p is plot- 
ted in Figure gj At threshold, emax/Tp [4(1 + 
m^/Arrip)]' 1 = 24.1%. At high energies, e max T p . 
Inverting this would give the minimum proton kinetic 
energy that makes a 7-ray photon with energy e, and 
restrict the integration in Equation ([24]) . Lacking a 
simple analytic inversion, it is simpler to limit the in- 
tegration numerically. 

The results of this part of study indicate that the 
resonance and scaling contributions to pion produc- 
tion are very different, at least for the function of 
Stephens & Badhwar [21] considered here. An im- 
provement to existing models could be obtained by 
adding different resonance contributions using the iso- 
bar model according to the cross sections shown in 
Figure [21 rather than transitioning between exclusive 
and inclusive formulations using empirical relations. 
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Figure 4: Kinematic maximum photon energy e max made in the process p + p — ► tt° — > 27, graphed as e macc /T p vs. T p . 



4. Determining the interstellar 
cosmic-ray spectrum 

Using the newly- measured emissivity spectrum [Bj], 
we can proceed to explore interstellar CR spectra that 
are compatible with it. The analysis is performed 
by first computing the matrices connecting model CR 
spectra to the observed gamma-ray emissivities in en- 
ergy bands, and then scanning the parameter space 
of the models. The method is Bayesian, allowing a 
complete scan of the parameters, computing poste- 
rior probability distributions, mean values and error 
bars correctly accounting for the correlations among 
the parameters. The method is Bayesian, allowing a 
complete scan of the parameters to compute poste- 
rior probability distributions, mean values, and error 
bars in order to correctly account for the correlations 
among the parameters. 

To illustrate the method, the measured emissivities 
have been fitted with 7-ray spectra calculated for a 
broken power law in momentum for protons and He- 
lium, with the free parameters being break momen- 
tum, spectral index below and above the break, and 
the overall normalization. The CR He/p ratio is fixed 
to the value measured by PAMELA at 100 GeV/nuc 
and the He and p spectra are assumed to have 
the same shape (since they cannot be distinguished 
in gamma rays). Note that CR spectra are expressed 
as particle density per momentum, for the reason ex- 
plained in [9]. The use of a sharp break in the CR 
spectrum is over-simplified but serves to illustrate the 
method; more physically plausible spectra are also be- 
ing investigated, but they do not lead to essentially 
different results. 

The first set of hadronic cross sections used is from 



Dermer [§[. The p-p cross sections are scaled for p- 
He, He-p and He-He interactions using the function 
given in Norbury & Townsend (TtJ. The second set 
is from Kamae et al. [l3| below 20 GeV, Kachelriefi 
& Ostapchenkojlll (QGSJET) above 20 GeV. For 
Kamae et al. [13] the p-p cross sections are scaled 
for p-He, He- p a nd He-He interactions using Norbury 
& Townsend \\J§. For KachelrieB & Ostapchenko [l2| 
the p-p, p-He, He-p cross-section are provided, so only 
He-He is scaled from p-p. The He fraction in the in- 
terstellar medium is taken as 0.1 by number. 

The electron (plus positron) spectrum producing 
bremsstrahlung is based on Fermi-LAT measurements 
above 10 GeV |2[, with a break below 3 GeV as indi- 
cated by synchrotron data [24[ . The synchrotron data 
shown there require a flattening of the interstellar elec- 
tron spectrum by about 1 unit in the spectral index 
below a few GeV, so this is used as a constraint; the 
actual low-energy index is determined by the gamma- 
ray fit. 

The resulting cosmic-ray proton and electron spec- 
tra and the corresponding emissivities are shown in 
Fig [5l The fit to the measured emissivities is good, 
as can be expected with the freedom allowed. It is 
encouraging that the proton spectrum is close to that 
measured directly at high energies, having been de- 
termined from gamma rays alone with no input from 
direct measurements except for the He/p ratio. The 
solar modulation is clearly seen in the deviation of the 
interstellar spectrum from the direct measurements 
below 10 GeV. Bremsstrahlung gives an essential con- 
tribution below ~ 1 GeV, and is an important com- 
ponent in the analysis. 

In this particular example, the interstellar proton 
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Figure 5: Spectra derived from model fitting. Emissivities from Casandjian 5], and cross sections from [(| (left 
column), [l2L (right column) . Model ranges are 1 standard deviation on the parameterized synthetic spectra. Top: 
Measured and derived cosmic-ray proton spectra. Data are AMS01 (asterisks) and PAMELA (diamonds). Yellow band 
shows model range. Middle: Measured and derived cosmic-ray electron spectra. Data are AMS01 (asterisks), PAMELA 
(diamonds), and Fermi-LAT (squares). Yellow band shows model range. Bottom: Fermi-LAT emissivity data (vertical 
bars) and model, with red and green curves showing the hadronic and leptonic bremsstrahlung contributions; the 
yellow band shows the total. 



spectrum 2 steepens by about 1/2 unit in the momen- 



2 Following the physical motivation explained in Q , we plot 
the density as a function of momentum n(p), because the flux 
j(T p ) = ^c/4tt)u(T p ) = ((3c/47r)\dp/dT p \n(p) = (c/4n)n{p). 



turn index above a few GeV, compatible with the ex- 
pectation from the cosmic-ray B/C ratio, which shows 
a similar break due to propagation. A power-law in- 
jection in momentum modified by propagation would 
then be a plausible scenario. The spectrum shown 
for Dermer [6| cross sections has momentum index 2.5 
(2.8) below (above) 6.5 GeV, with a scaling factor 1.4 
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relative to PAMELA at 100 GeV; for the Kamae et 
al. [l3|, KachelrieB & Ostapchenko [Hj cross sections 
the values are 2.4 (2.9) and 1.3, with the same break 
energy. The sensitivity of the results to the cross sec- 
tions is evidently significant but not overwhelming, 
though the range of uncertainties in the production 
cross sections must be included in the error budget 
to derive firm conclusions. In both of these illustra- 
tive cases, the high-energy proton index is compatible 
with PAMELA (2.82) [4]. The scaling factor excess 
may have various origins, including uncertainties in 
the cross sections and the gas tracers, hidden system- 
atic errors in the direct measurements themselves; a 
combination of these is possible. A difference between 
the interstellar spectrum and the direct measurements 
cannot be ruled out at this stage either. 

This is an example of how it will be possible to 
constrain the interstellar CR spectra with the Fermi- 
LAT emissivity data. In a forthcoming paper [25|, all 
the uncertainties will be addressed, including those 
in emissivities (gas, instrumental response, etc.) and 
cross sections. The evidence for a break and more 
exact constraints on the spectrum will be obtained 
there. 

The authors wish to thank Michael Kachelriess and 
Sergey Ostapchenko for valuable discussions on the 
use of their cross-section code. 
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